Exactly Integrable Dynamics of Interface between Ideal Fluid and Light Viscous Fluid 



Pavel M. Lushnikov 1 ' 2 

1 Theoretical Division, Los Alamos National Laboratory, MS-B284, Los Alamos, New Mexico, 87545 
2 Landau Institute for Theoretical Physics, 2 Kosygin Str., Moscow, 119334, Russia^} 

It is shown that dynamics of the interface between ideal fluid and light viscous fluid is exactly inte- 
grable in the approximation of small surface slopes for two-dimensional flow. Stokes flow of viscous 
fluid provides a relation between normal velocity and pressure at interface. Surface elevation and 
velocity potential of ideal fluid are determined from two complex Burgers equations corresponding 
to analytical continuation of velocity potential at the interface into upper and lower complex half 
planes, respectively. The interface loses its smoothness if complex singularities (poles) reach the 
interface. 
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Dynamics of an interface between two incompressible 
fluids is an important fundamental problem which has 
numerous applications ranging from interaction between 
see and atmosphere to flow through porous media and 
superfluids. If one neglects gravity and surface tension, 
that problem can be effectively solved in some partic- 
ular cases in two dimensions with the use of complex 
variables. Integrable cases include Stokes flow of viscous 
fluid with free surface , dynamics of free surface of ideal 
fluid with infinite depth and finite depth 0, dynam- 
ics of an interface between two ideal fluids , ideal fluid 
pushed through viscous fluid in a narrowgap between 
two parallel plates (Hele-Shaw flow) [a, la 13 • 

Here a new integrable case is found which corresponds 
to two-dimensional motion of the interface between heavy 
ideal fluid and light viscous fluid in absence of gravity 
and capillary forces. The interface position is given by 
z = rj{x,t), where the first, heavier fluid (indicated by 
index 1) with the density p\ occupies the region — oo < 
z < rj(x,t) and the second, lighter fluid (index 2) with 
the density p 2 occupies the region r)(x,t) < z < oo. 

Suppose that the kinematic viscosity of the fluid 2, u 2 , 
is very large so that fluid's 2 flow has small Reynolds 
numbers and, neglecting inertial effect in the Navier- 
Stokes Eq., one arrives to the Stokes flow Eq. 

z/ 2 V 2 v 2 - — Vp 2 = 0, V • v 2 = 0, (1) 

P2 

where v 2 = (t> 2 . x , w 2 , z ) is the velocity of the fluid 2, V = 
(d x ,d z ), and p 2 is the fluid's 2 pressure (similar physical 
quantities for the fluid 1 have index 1 below) . Additional 
assumption necessary for applicability of Eq. is a 
small density ratio, 



p 2 /pi<l, pi = l, 



(2) 



which ensure that the fluid 2 responds very fast to per- 
turbations of the interface as inertia of the fluid 2 is very 
small compare with fluid's 1 inertia while time dependent 
perturbations of the fluid 2 decay very fast due to large 
viscosity v 2 . According to Eq. (JTJ, the response of the 



fluid 2 to motion of the interface is static. For any given 
normal velocity of the interface, v n , Eq. allows to 
determine the pressure p 2 \ z=v at the interface. In other 
words, the fluid 2 adiabatically follows the slow motion 
of the heavy fluid 1 and Reynolds number of the fluid 2 
remains small at all time. 

The velocity of the potential motion of ideal the fluid 1, 
vi = V</), can be found from solution of the Laplace Eq., 
V 2 = 0, which is a consequence of the incompressibility 
condition, V • Vx = 0, for potential flow. Boundary con- 
ditions at infinity are decaying, |vi|, p\ — > for z — > 
— oo; 

|v 2 |,p 2 — >0 for z — > +oo. 

Motion of the interface is determined from the kine- 
matic boundary condition of continuity of normal com- 
ponent of fluid velocity across the interface: 



v 2n \ z = n =dtv[l + (d x r)) 2 ] 1/2 , (3) 
and n= (-9,77, mi+ta^r 172 



where v 1{2 )n = n-v 1(2) 
is the interface normal vector. 

A dynamic boundary condition is a continuity of stress 
tensor, a 1{2) , jm = -p 1{2) 6 jm + <y' 1[2) ]m , <j' m jm = 

Pi(2)^i(2)(— b 9x ^ ), xi = x, x 2 = z, across 

the interface: rijUi j m \ z=v — nja 2 , j m \ z=v (repetition of 
indexes j, m means summation from 1 to 2), which gives 
two scalar dynamic boundary conditions: 



Pi\z= n =P2\z= n +n m nj(r 2tjm \ z =n, 
lmnjO- 2 j m \ z=rj = 0, 



(4) 



where the absence of viscous stress in the ideal fluid 1, 
v\ = 0, is used, n m ,l m are components of the interface 
normal vector, n, and the interface tangential vector, 1 = 

(l,d x r])[l + {d x rj) 2 ] 172 . The pressure p\ of the fluid 1 
at the interface can be determined from a nonstationary 

Bernoulli Eq., [d t 4> + ±(V0) 2 + f] =0. 

To obtain a closed expression for interface dynamics 
in terms of fluid's 1 variables only, one can first find an 
expression for the pressure at the interface through the 
normal velocity v n . 
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It follows from Eq. JIJ that V 2 p 2 — and the 
Fourier transform over x allows to write the solution 
of the Laplace Eq. with the decaying boundary con- 
dition at x — > oo as P2k( z ) — P2fc(0) exp(— \k\z) = 
J dxp2{x,z) e~xjp(—ikx). 

To determine ~V2\z=r) one can introduce a shift oper- 
ator, L 2 , defined from series expansion: v 2 (x, z)\ _ = 



L 2 Mx,0) = [(l+vdz + ^rfd 2 z + . . .)v 2 (ar, z) 
Eq. Q to find w 2 , xfc (z) = [c fe - ikz^gfa 
V2,zk{z) = isign(fc)cfe + (|/c|z + 1) 



P2fc(0) 



and use 

exp(-\k\z), 

<2 P2V2 \k\\ exp(-|/c|z), 
where V2,x k (z), V2,zk{z) are the Fourier transform over 
x of the components of the velocity v 2 and functions Ck, 
P2 k (0) should be determined from the dynamic boundary 
conditions JIJ. 

Operator X 2 can be expressed, using Eq. l|T)l. in terms 
of the operator fe : L% = 1 — ?7fc + \rfk + . . . , where the 
integral operator k is an inverse Fourier transform of \k\ 
and is given by 



ax 



(•5) 



i 2 9 2; v 2 (a;,0), [9 z u 2 ,x(a;, z)]L _ r 
L 2 fc" 1 9 :!; p 2 (a;,0), p 2 {x,r]) 



2f2 P2 



Here #7(x) = ±P.V. f + °° 4^cte' is the Hilbert trans- 

•> \ / 7T J — OO X ' — X 

form and P.V. means Cauchy principal value of inte- 
gral. H can be also interpreted as a Fourier transform of 
i sign(fc). 

In a similar way one can show that [d x V2(x, z)]\ z _ = 

= -Lkv 2 , x (x,0) - 
L 2 p 2 (a;,0) and using 
kinematic ((3J) and dynamic Q boundary conditions one 
can find Pi(x, rj) as a linear functional of v n . That lin- 
ear functional can be expressed in a form of powers series 
with respect to small parameter | c?^ 77 1 , which has a mean- 
ing of typical slope of the interface inclination relative to 
the interface undisturbed (plane) position. 

At leading order approximation over small parameter 
\d x T)\ one gets: pi(x,r}) = p2(x,0), v n = V2,z(x,r]), and, 
respectively, response of pressure to normal velocity is 
given by 



Pi I 



2p 2 V2kv„ 



(6) 



In other words, Eq. © determines a static response of 
the fluid 2 to the motion of the interface. 

Eq. JfJJl together with the kinematic boundary condi- 
tion © and the Laplace Eq. for the velocity ponetial <j> 
completely defines the potential motion of the fluid 1. 

Following Zakharov |£g, one can introduce the surface 
variable tp(x) = 4>{x, 77), which is the value of the velocity 
potential, cj)(x,z), at the interface. Kinematic boundary 
condition @ can be written at leading order over small 
parameter \d x r]\ as 



where a new function, v — d x ip, is introduced which has 
a meaning of the tangent velocity of the fluid 1 at the 
interface. 

Similar to the shift operator L 2 , one can define a shift 
operator, L\ = l + rjk+^^k 2 + . . . , which corresponds to 
the harmonic function </> with vanishing boundary condi- 
tion </> — > for z — > — 00. A Fourier transform of 0, 
4>k(z) — (f>k{0) exp(|fc|z), allows to find the components of 
fluid velocity at the interface: (d x (j))\ _ — Lid x <j)(x, 0) = 

I Z — 1] 

LidxLT 1 ^, (d z (f>)\ _ — LikLT 1 ^ through surface vari- 

I z — rj 

ables T], ip. Time derivative (j>t in the nonstationary 
Bernoulli can be found from d t ip = d t 4>\ z =ri + dtV^z^lz^ 
and one gets at leading order approximation over | c?^ 77 1 : 



dv 
di 



-A<: 



(Hv) 



2v2P2d 2 x v. 



(8) 



Note that Eq. JHJ does not include variable rj which is a 
peculiar property of lowest perturbation order over Ic^ryl- 

Because the surface tension and gravity is neglected 
here, the total energy of two fluid equals to total kinetic 
energy, K. K decays, ^ ~ f dfc|fc| 3 |^| 2 < 0, due 

to dissipation in the fluid 2. If the fluid 2 is absent, which 
corresponds to p 2 = 0, then K is conserved, =£■ — 0, and 
the motion of the fluid 1 can be expressed in the standard 
Hamiltonian form 0,0: §? = , §£ = 

Equations, similar to (JJJ) , (jHJ can be derived for three 
dimensional motion also with the main difference that the 
operator k in three dimensions is not given by (JSJ) but 
determined from the Fourier transform of |fc| over two 
horizontal coordinates. Subsequent analysis is however 
restricted to two dimensional fluid motion only. 

The real function v(x) can be uniquely represented 
as a sum of two complex functions ?/ + ) and v^~\ v — 
[«(+) +t;(~)]/2, which can be analytically continued from 
real axis x into upper and lower complex half-planes, re- 
spectively. The Hilbert transform acts on these functions 
as Hv {+) = iv {+ \ Hv { - ] = -iv^ and Eq. © splits 
into two decoupled complex Burgers Eqs. for and 
„(-)• 



at 



(9) 



d t r) = -Hv, 



(7) 



where an effective viscosity, v = 2^ 2 p 2 is introduced 
to make connection with the standard definition of real 
Burgers Eq. H3]. Similar reduction of integro-differential 
Eq. (like Eq. (JSJ) to complex Burgers Eq. was done in 
Ref. 

If the fluid 2 is absent, v — 0, complex Burgers Eqs. 
@ are reduced to inviscid Burgers Eqs. (the Hopf Eqs.) 
which were derived for ideal fluid with free surface in Ref. 
(note that definition of in this Letter differs from 
similar definition in Ref. @ by a factor 1/2). While 
viscosity V2 is large enough to make sure that Reynolds 
number in the fluid 2, R 2 , is small, R 2 ~ v n /{v2k) <C 1 
(k is a typical wave vector of surface perturbation) but 
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effective viscosity v can be smaii provided p2 <C R2 <C 1 
so that Reynolds number , R, in complex Burgers Eq. 
@ is large, R ~ R 2 /p 2 > 1. 

Complex Burgers Eq. is transformed into the com- 
plex heat Eq. 9w g t - = vd^u^ via the Cole-Hopf trans- 
form -y (±) = -29^^-. Solution of the heat Eq. 
with initial data (x, t) | . = u^\x), u^(x,t) — 
{A-Kvt) -1 / 2 dx' exp — ^ Uq ± '(x / ), is an ana- 



lytic function in complex x plane for any t > because 
integral of right hand side (rhs) of this Eq. over any 
closed contour in complex x plane is zero (Morera's the- 
orem). Then, according to the Cole-Hopf transform, so- 
lution of the complex Burgers Eq. can have pole sin- 
gularities corresponding to zeros of iv^'(x,t). Number 
of zeros, n(j), of u^'(x,t) (each zero is calculated ac- 
cording to its order) inside any simple closed contour 7 
equals to / dxd x u^' (x , t) / u^' (x , t) . Integration of 
Eq. © over 7 allows to conclude that 71(7) is conserved 
as a function of time provided zeros do not cross 7. Thus 
number of zeros in entire complex plane can only change 
in time because zero can be created or annihilated at 
complex infinity, x — 00, provided u^'(x, t) has an es- 
sential singularity at complex infinity. 

From physical point of view it is important that zeros 
oiu^(x,t) can reach real axis x — Re(x) which distin- 
guishes the complex Burgers Eq. from the real Burgers 
Eq. Solution of the real Burgers Eq., which corresponds 
to Eq. © with v^(x,t)\ t=Q = Re[v {± \x, t)\ t=Q \, has 
global existence (remains smooth for any time), while so- 
lution of the complex Burgers generally exists until some 
zero of u^\x, t) hits real axis x for the first time. 

To make connection with inviscid case one can look 
at initial condition for (x, 0) with one simple pole in 
the lower half-plane: 



.(+) 



Or,0) 



2A 

x + ia 



Re(a) > 0. 



(10) 



Solution of the inviscid {p = 
initial condition (|10f> gives 



0) Burgers Eq. with 



4.4 



2]: 



v 



(+) 



inviscii 



,(x,t) 



— , which has two moving branch points: 

x+ia+V (x+ia) 2 -SAt 

Xi t 2 = —ia± 2y '2 At. One of these branch points reaches 
real axis in a finite time if either A < or Re(A) 7^ 0. 
As the branch point touches the real axis, the inviscid 
solution is not unique any more and the interface looses 
its smoothness 0- 

Consider now solution of the viscous Burgers Eq. © 
with nonzero effective viscosity v and with the simple 
pole conditions 1)10(1 . Respectively, initial condition for 



the heat Eq., is given by u Q 



(+) 



x + ia) A ' v and has 
branch point at x = — ia. Solution of the heat Eq. 

gives w^ + ^(x,i) = exp (i^)if^^ — ^^[x + ia]\, where 

p, = — A/v and H^{z) is the Hermite function defined as 



f °°dye-yy» 



cos 



(2zy 



Zeros of 




FIG. 1: Motion of poles of the complex velocity, u' + '(:E,t), in 
complex x plane for A = —1/8, ft = 7.6, v ~ 0.0164, a = 
1, t = 0.8. Arrows point out the position, direction and 
magnitude of moving poles. Uppermost arrow designates 
the pole which corresponds to zero of the Hermite function 
with the largest real part (that pole first reaches real axis 
for t V i SCOUS ~ 1.91 producing singularity of the interface sur- 
face). Dotted line connects two branch points (filled circles) 
of inviscid solution. Upper branch point reaches real axis for 
tinvisdd = 1 which corresponds to singularity in the solution 
of the inviscid Burgers Eq. Empty circle designates the simple 
pole initial condition 11101 . Viscous solution becomes singular 
at later time compare to inviscid solution, t V i 3COUS > tinvisdd- 



v,( + '(x,t) (and, equivalently, poles of v^ + \x,t)) move in 
complex x plane with time as (see Figure 1) 



Xj (t) = i{2y vtzj — a) 



(11) 



where z\, Z2, ■ ■ ■ , are complex zeros of the Hermite func- 
tion. 

Consider a particular case, jl = n, n is a positive inte- 
ger number. The Hermite function is reduced to the Her- 
mite polynomial H n {z) which has n zeros, Z\, Z2, ■ ■ ■ , z n 
located at real axis z = Re(z), z n corresponds to the 
largest zero. Location of real zeros of the Hermite func- 
tion with real jl is close to location of zeros of the Hermite 
polinomial with the closest integer n to the given jl while 
zeros with nonzero imaginary part (which corresponds 
to tails with nonzero real part in Figure 1) disappear for 
jl = n. Zeros of the Hermite polinomial are moving with 
time parallel to imaginary axis x = Im(x) in complex x 
plane according to l|ll|) and the complex velocity is 
described by set of moving poles: 



,(+) 



i 



Xj(t)' 



(12) 



v^' is given by the same expression with Xj replaced by 
their conjugated values Xj. 

Eqs. (O have also another wide class of solutions, 
"pole decomposition", corresponding to Eq. (|12f) with 
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X 

FIG. 2: The interface position, rj{x, t), according to solution 
of Eqs. J7J, with finite viscosity, v — 1/64 (solid line) and 
zero viscosity, v = 0, (dotted line) for A — —1/8, a — 1, t — 
1. Viscous solution has 8 moving poles while inviscid solution 
is singular at x = (d£ x T}\x=o — > — oo as £ — » t inviaM = 1). 
Both solutions are almost indistinguishable outside a small 
neighborhood around x = 0. As /x increases, the viscous so- 
lution approaches invisicid. 



= — 2v — z — i 71 is arbitrary positive integer 

|l3|. Simple pole initial condition (|10|) with /i = n is 
particular case for which Xj\ t =o = for any j. 

As f(a;,t) is known from solution of the heat Eq. and 
the Cole-Hopf transform one can find rj{x, t) from Eq. 
(|7|). Interface dynamics is determined from the most 
rapid pole of v^' which first reaches real axis, x = Re(x). 
E.g., for initial condition l)10[l. the pole singularity of 
first hits real axis, x = Re(x), from below at time 

tviscous = 40Re{zl a:c y^ wnere z max is a complex zero 
of the Hermite function with the largest real part for 
given jl. Simultaneously, the pole singularity of first 
hits real axis from above at the same point. Figure 2 
shows rj(x) at the time, t — ti nv i SC id, when singularity 
(branch point) of inviscid solution first reaches the inter- 
face breaking analyticity of inviscid solution. It is seen 
that viscous solution significantly deviates from inviscid 
one only in the narrow domain around x = 0. 

Viscous solution remains analytic for t > ti nv i SC id un- 
til t < t viscous (t viscous ~ 1.91 for parameters in Fig.l). 
However, for t — > t V i scous , surface elevation behaves as 
f] ~ (-a/2z% lax ) log [x 2 + (2VWz max - a) 2 ] near x = 
(it is set here Im(a) = Im(A) — 0) meaning that small 
slope approximation used for derivation of Eqs. 0, (JSJ) 
is violated for t — > t V i SCOUS and full hydrodynamic Eqs. 
should be solved near singularity. One can find a range 
of applicability of Eqs. (01, © by looking at correction 
to these Eqs. E.g. the analysis for parameters of Fig. 1 
shows that the correction is important for t > 0.9t V i SCOUS 
(for t = 0.9t V i SCOUS correction to 77| x =o is about 30%.) 
Detail consideration of that question is outside the scope 
of this Letter. Note that the question whether an actual 
singularity of the interface surface occurs in full hydro- 
dynamic Eqs. remains open. 



To make connection with dynamics of ideal fluid with 
free surface (corresponds to the inviscid Burgers Eqs.) 
one can consider a limit v — > and, respectively, jl — > oo. 
It can be shown from the asymptotic analysis of the in- 
tegral representation of the Hermite function that the 
largest zero, z max , is given by z max = 2 1 / 2 /i 1 / 2 +0(/i~ 1 / 6 ) 
The leading order term, 2 1//2 /i 1 / 2 , exactly corresponds to 
the position of the upper branch point of inviscid solu- 
tion (see Fig. 1) while term 0(/i -1 / 6 ) < is responsible 
for the difference between ti nv i SC i,i and t V i SCOUS . Even for 
moderately small v as in Fig. 1 that difference is numer- 
ically close to 1 because of small power pT 1 ^ . 

It is easy to derive a wide class of initial conditions 
for which solution Q , (JBJ exists globally and the inter- 
face remains smooth at all times. E.g. one can take 
= aoe^^^o 1 , fc = Re(k ) > or any sum of 
imaginary exponent which ensure that there is no zeros 
at Im(x) = 0. However, it we suppose that there is a 
random force pumping of energy into system (or random 
initial condition) then one can expect that some trajecto- 
ries with nonzero measure would have poles which reach 
real axis in a finite time. 

In conclusion, one can mention possible physical ap- 
plications. Eqs. JTJl, © describe a free surface dynamics 
of Helium II with both normal (yi ^ 0) and superfluid 
(y\ = 0) components. Derivation of these Eqs. is slightly 
different from given in this Letter because both fluids oc- 
cupy the same volume but resulting Eqs. are exactly the 
same as 0), JSJ). For classical fluids viscosity is nonzero 
but v\ can be neglected and the fluid 1 can be considered 
as ideal fluid provided the ratio of dynamic viscosities of 
two fluids is large, vipij(y\p\) ^> 1. E.g. that ratio is 
~ 900 for glycerin and mercury while ratio of their den- 
sities is ~ 0.09 which makes them good candidates for 
experimental test of the analytical result of this Letter. 
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